function [CL,CD] = getCLCD(alpha, Ma)
    alphadeg = alpha*180/pi;
	p00 = 0.5956;  p10 = -0.09572; p01 = -0.06861; p20 = 0.0005887;
	p11 = 0.01532; p02 = 0.003406; p30 = 0.0005817; p21 = -0.00101;
	p12 = -0.000579; p40 = -3.246e-05; p31 = 2.228e-05; p22 = 3.615e-05;
	p50 = 5.154e-07; p41 = -5.753e-08; p32 = -7.175e-07;
	x = Ma;
	y = alphadeg;
	CD = p00 + p10*x + p01*y + p20.*x.*x + p11.*x.*y + p02.*y.*y + p30.*x.*x.*x + p21.*x.*x.*y...
		 + p12.*x.*y.*y + p40.*x.*x.*x.*x + p31.*x.*x.*x.*y + p22.*x.*x.*y.*y + p50.*x.*x.*x.*x.*x...
         + p41.*x.*x.*x.*x.*y + p32.*x.*x.*x.*y.*y;
    p00 = -0.04928;		p10 = -0.009818;	p01 = 0.05259;		p20 = 0.0011;
    p11 = -0.001381;	p02 = 0.0004019;	p30 = -2.586e-05;	p21 = 2.176e-05;
    p12 = -2.289e-06;
	CL = p00 + p10.*x + p01.*y + p20.*x.*x + p11.*x.*y + p02.*y.*y + p30.*x.*x.*x +...
         p21.*x.*x.*y + p12.*x.*y.*y;
end